Topological interpretation of subharmonic mode locking in coupled oscillators with 

inertia 
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A topological argument is constructed and applied to explain subharmonic mode locking in a 
system of coupled oscillators with inertia. Via a series of transformations, the system is shown to be 
described by a classical XY model with periodic bond angles, which is in turn mapped onto a tight- 
binding particle in a periodic gauge field. It is then revealed that subharmonic quantization of the 
average phase velocity follows as a manifestation of topological invariance. Ubiquity of multistability 
and associated hysteresis are also pointed out. 
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In a nonlinear oscillator system driven periodically, the 
competition between the natural frequency and the driv- 
ing frequency in general leads to either an almost periodic 
motion or a periodic one, depending on the parameter 
range Q|. The latter, called mode locking, is character- 
ized by the quantization at rational values of the average 
phase velocity. In particular subharmonic mode locking, 
appearing in the presence of the inertia term, results in 
the devil's staircase structure. One of the well-known ex- 
amples is the Josephson junction driven by combined di- 
rect and alternating currents, with the capacitance play- 
ing the role of inertia . Governed by the same equation 
of motion as a driven pendulum, it displays dc voltage 
plateaus in the current- voltage characteristics, known as 
Shapiro steps jj). Similar voltage quantization has also 
been observed in arrays of Josephson junctions, yield- 
ing integer giant Shapiro steps and subharmonic 
steps H according to the absence/presence of capacitive 
terms. Unfortunately, in spite of the deceptively sim- 
ple equation of motion, even the smgfe-junction problem 
has resisted complete analytical solutions, especially, in 
the presence of the capacitive term, except for the re- 
sults mainly based on the circle map j?]] and on the ap- 
proximate analysis by means of expansion and averag- 
ing m. Accordingly, such mode locking phenomena in 
arrays have been demonstrated mostly by numerical sim- 
ulations. On the other hand, the topological argument, 
proposed for the system without the capacitive term [fl0[ , 
reveals topological invariance of the system as the phys- 
ical origin of quantization p| . As in the case of the 
quantum Hall effect |T^ ], the topological argument does 
not provide quantitative information, e.g., on the locking 
structure. Nevertheless it not only clarifies the nature of 
quantization but also provides a link between dynamics 
and statics by interpreting (dynamical) mode locking in 
terms of (static) topological invariance. 

In this work, we construct a topological argument for 
the system with inertia, and apply the idea to Josephson- 
j unction arrays or systems of coupled oscillators, with 



attention to the resulting subharmonic locking. For this 
purpose, we consider an appropriate canonical transfor- 
mation of the dynamic equations of motion and the cor- 
responding Fokkcr-Planck equation, the stationary solu- 
tion of which gives the effective Hamiltonian in the form 
of a classical XY model with periodic bond angles. Via 
mapping onto a tight-binding particle in a periodic gauge 
field, subharmonic quantization of the average (dc) phase 
velocity is revealed as a manifestation of the topological 
invariance. Also suggested is ubiquity of multistability, 
providing a natural explanation of the observed hystere- 
sis due to the inertia. 

We begin with the set of equations of motion for a 
system of N coupled oscillators 



N 



^2[Mij(f>j + jMijtpj + Jij sin(fa - <j>j)} = h 



(1) 



where fa represents the phase of the ith oscillator, M^- 
the (rotational) inertia matrix, 7 the damping parameter, 
and measures the coupling strength between the oscil- 
lators i and j. The right-hand side describes the periodic 
driving with frequency Q: U = Ii t d+Ii,a cos fii, where Ii t a 
and Ii :a are the amplitudes of the dc and ac components, 
respectively, of the driving on the ithe oscillator. There 
are two cases depending on the detailed form of the iner- 
tia matrix. The simple case that My = M£y describes 
the system of coupled oscillators, each of which possesses 
inertia M and suffers from dissipation of strength 7M un- 
der driving Ii. On the other hand, with fa denoting the 
phase of the superconducting order parameter at site i, 
Eq. (0) describes the dynamics of the array of resistively 
and capacitively shunted junctions (RCSJs), where the 
combined direct and alternating current Ii is fed into the 
grain at site i ||[l3|. In this case, My corresponds to the 
capacitance matrix and assumes the form My = CAy 
with the junction capacitance C and the lattice Laplacian 
Ay = zSij—Sij/ , where j' represents the neighboring sites 
of j and z the number of such neighbors. The damping 
parameter 7 is inversely proportional to the shunt resis- 



1 



tance of the junction. For simplicity, we henceforth con- 
centrate on the case My = M<5y since the generalization 
to the case of My = CAy is straightforward. 

Equation (Q) may be written in the form of Hamilton's 
canonical equations: fa = dH/dpi and & = —dH/d(f>i 



with the Caldirola-Kanai Hamiltonian 14 



h C k = ^ E^ + Qi) 2 ^ 1 - E J ^ cos (^ - ^) e7 *' 



(2) 



where (fa,Pi) are conjugate variables and Qi is the 
"gauge charge" given by Qie~ 7t = U or 



7 



7 2 + fi : 



r e 7t (7COsrtt + ftsinfit) + Q° (3) 



with arbitrary constant In this sense the classical me- 
chanical system with dissipation, governed by the equa- 
tions of motion can be described by the Hamilto- 
nian in Eq. (0) flq| . We then introduce new variables 
(9i,Pi) according to pi = d^/dfa and 6>; = d<&/dpi with 
the generating function $({pi},{fa}) = ^2iPi{fa + a i), 
where Men = —e ft Qi. Under this canonical transfor- 
mation, Eq. (^|) yields, apart from a constant term, the 
new Hamiltonian 



Hck = Hck + -^r 



d<S> 
~dt 



i<j 



Jije 11 cos(6i - 9j - ay), (4) 



where pi (= Pi) is conjugate to In view of the RCSJ 
array, where a uniform current is usually fed into the 
sites along one edge and extracted from those along the 
opposite edge, we consider the case that some oscillators 
are driven by / = Id + I a cos fit and some others by — /. 



Then the bond angle ay = 
either vanishes or becomes 



dj, depending on 



= T 



M 



7 



-t + 



7 



Ja n 

2 + ft 2 \fi 



sin fit — cos fit 



(5) 



apart from an arbitrary constant. 

Note that the energy of the system is given by 
e -i t H CK p5| ; this also corresponds to the effective 
Hamiltonian describing the statistical mechanics of the 
system. To see this, we for the moment consider the sys- 
tem at finite temperatures and generalize the equations 
of motion (|l|) appropriately: 



N 



^jM^fa + jM lj fa + 

.7=1 



Jij sin( 



fa)]=I i + m , (6) 



where r\i is the random (thermal) noise acting on the ith 
oscillator. In the system of oscillators with My = M£y, 



the noise is characterized by the zero mean and the cor- 
relations 

= 2 1 Mk B T5 l3 5{t-t l ) 

at temperature T. In the RCSJ array with My = CAy, 
f]i is given by the sum of the noise currents from neigh- 
boring sites: rji = V ; '/,/, ,■ with 

(%(*)%;(0> = 2 1 Ck B T(5 lk d, l ~ 5u5 jk )5(t-t'). 

Motivated by the canonical transformation in the ab- 
sence of noise, we write Eq. (0) in the form 



= M-V*e" 7 * 



Pi 



e 7 ' sin( 



ay) + ?7ie 7 *, 



(7) 



where 9i = fa + 0.$. The set of Langevin equations (J? 
may be transformed into the Fokker-Planck equation |l1 



dP 



dt ~ 2s M 6 d9i 



Jy-e 7 * sin( 



h — 9j — aij)-^— 



dP 

dpi 



-fMkBT^e 21 



dpi ' 



(8) 



which describes the time evolution of the probability dis- 
tribution P ({#;}, {pi}, t) of phases and momenta at time 
t. Equation (g) yields the stationary solution valid in the 
limit t — > oo: 

P({^},fe})(x e -^'/ fe « T , 
where the effective Hamiltonian is given by 



H eff ~ 



-^^Tpf - y " cos(9i 



i<j 



(9) 



with M = Me 2lt . It is thus concluded that the sta- 
tionary distribution has the form of a Gibbs measure, 
with the effective Hamiltonian indeed corresponding to 
the energy of the system {Heff — e~ Jt HcK)- 

The first term in Eq. (M becomes vanishingly small in 
the stationary state (t — > oo); it is further obvious that 
the kinetic energy in the above classical system decouples 
from the interaction energy. We thus obtain the classical 
XY Hamiltonian 111 



H 



= - J E 



cos 



- 9j - ay), 



(10) 



where the nearest-neighbor coupling (Jy = J5y) has 
been assumed for convenience. At zero temperature, 
which is our concern, the system described by the Hamil- 
tonian ( |io| ) is equivalent to a tight-binding particle (of 
charge e), with 2ksT/J taking the role of the energy 
eigenvalue jTJj) . The Hamiltonian describing such a tight- 
binding system has the position representation 



(i\H\j)=e 



—laa r 



(11) 
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where |z) is the position eigenket and ay may be viewed 
as the line integral of the appropriate gauge potential a: 
a ij = ( e /^ c ) // el cII. Equation (^) shows that ay, defined 
modulo 2tt, is periodic in time with period r = 27rm/f2 
only if 



i d = -jMn 

to 



(12) 



with to and s integers. 

When Eq. ( p"2| ) is satisfied, the Hamiltonian ( [Tl| ) as 
well as the gauge field ay has the periodicity r and 
the Floquet theorem is applicable to the corresponding 
Schrodinger equation for the wave function ^ = (i| v E'), 
giving the relation pOfl 



-iEt 



(13) 



with the quasienergy E. This imposes that, apart 
from the dynamical contribution Et, the corresponding 
change in the phase of the wave function \&j should be 
an integer multiple of 2ir: 9i(t+r) — 6i(t) = 2riiir — Et, 
where n, is an integer depending on the form of ay, i.e., 
of the driving li and E has been assumed to be real (see 
below). We thus have the average change rate of the 
phase 



I r n- 
i) = - dt9i= —n-E 
T Jo m 



or in terms of the original phase 



I-Ld 

7M 



E, 



(14) 



(15) 



where hi = m±s or m for driven (Ii y d = Hd) or undriven 
(Ii,d — 0) oscillators, respectively. Accordingly, the av- 
erage change rate of the phase difference or the average 
(dc) phase velocity, which usually gives the appropriate 
physical quantity, e.g., the voltage in the case of an RCSJ 
array, indeed displays subharmonic mode locking: 



(Vy) = <&) 



n „ 

— n 

TO 



(16) 



with n = hi — hj. Note that for given configuration of 
driving the integer n in Eq. ( |l6| ) is determined by the 
winding number rij, manifesting the topological nature 
of the mode locking. 

The subharmonic mode locking given by Eq. ( [l6| ) can 
persist even for the (dc) driving slightly off the condition 
in Eq. (|2h. To see this, we take 



Id 



s 

m 



(17) 



for small e, which leads to the bond angle either zero or 
fly = a®j =F eQ,t with a^- representing the periodic part 
for e = 0. In this case the Hamiltonian of the system as 
well as ay is in general not periodic. However, in terms 



of the shifted phase Xi = @i + e»Ot, where ej = ±e or 
for driven or undriven oscillators, the Hamiltonian in 
Eq. @ reads 



H 



XY 



-j c ° s (xi 



Xj 



(18) 



where periodicity has been restored. Accordingly, the ar- 
gument leading to Eq. ( |l4| ) is applicable, giving (xi) = 
(rii/m)Q — E. This in turn leads to 



i) = 0i) + 



—n 

TO 



H,d 

E, 



(Xi 



EiQ, ■ 



1 i,d 

jM 



(19) 



which is precisely Eq. (|l5|). It is thus obvious that 
Eq. (16) remains unchanged: Given subharmonic quan- 
tization can persist in the finite interval of the dc driving 
strength Sid oc e, thus generating the step structure in the 
appropriate response characteristics. Since rational num- 
bers form a dense set, this also indicates that there can 
exist multiple quantization states with different values of 
n/m in Eq. ( |l6| ) for given driving strength Id- However, 
some of those, corresponding to complex values of the 
quasienergy E with nonzero imaginary parts, are unsta- 
ble J|,^(|, and only stable states (with real values of E) 
among those determine the step structure in the actual 
response characteristics. Although such information as 
stability cannot be obtained by the topological argument, 
it is certainly plausible to have two or more states stable 
in some intervals of the driving strength. Such multista- 
bility in general gives rise to hysteresis behavior in the re- 
sponse characteristics, which has widely been observed in 
the oscillator systems with inertia terms |^^| JT7| , pl[ . Note 
also that the topological argument is in essence a zero- 
temperature analysis like in Refs. |l^,|l^]: At finite tem- 
peratures thermal fluctuations make the mapping onto 
Eq. (|ll]) and the following argument inexact. Neverthe- 
less, at sufficiently low temperatures, the phase slippage 
induced by fluctuations should have an exponentially low 
rate, hardly affecting the quantization itself determined 
by the winding number. On the other hand, it is expected 
that fluctuations tend to destabilize various quantization 
states, reducing the stability interval and the emergence 
of multistability. This is consistent with smoothing out 
of the step structure Q| and suppression of hysteresis p2| , 
observed in the presence of noise. 

In summary, we have constructed a topological argu- 
ment for subharmonic mode locking in a driven system 
of coupled oscillators with inertia. Starting from the dy- 
namic equations of motion, we have derived the effec- 
tive Hamiltonian for the system in the form of a classical 
XY model with periodic bond angles, which in turn has 
been mapped onto a tight-binding particle in a periodic 
gauge field. It has then been shown that subharmonic 
quantization of the average (dc) phase velocity follows as 
a manifestation of topological invariance, revealing the 
topological nature of the dynamical mode locking. Also 
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revealed is the possibility of multistability, providing a 
natural explanation of the ubiquity of hysteresis due to 
the inertia. In view of the fact that the set of equations 
of motion (1) describes a prototype of oscillatory system, 
we believe the result of this paper to be rather general 
and applicable to a variety of oscillatory systems display- 
ing mode locking phenomena. 
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